Coherent phase shps in superconducting nanorings 
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We study quantum fluctuations of persistent current in a small superconducting ring. Based on a 
microscopic model of the ring we argue that under certain conditions such ring will exhibit coherent 
quantum phase slips, similar to those in a flux qubit. We evaluate the frequency of such coherent 
oscillations and flnd that it is strongly dependent on wire's diameter primarily due to the large 
momentum released by the condensate as a result of a phase slip event. We also find that the value 
of such frequency is not a self-averaging quantity, that is, it depends on a particular realization of 
the static impurity potential. 

PACS numbers: 74.78.Na, 74.40.-n, 74.50.+r, 85.25.-j 



Relaxation of persistent currents in quasi one- 
dimensional (ID) superfluids due to phase slip transi- 
tions has been a subject of extensive studies over several 
decades. While at finite temperatures phase slips are 
known to occur via thermal activation through an ef- 
fective barrier [l|, at sufficiently low temperatures they 
are believed to be generated by the quantum fluctua- 
tions of the order parameter, i.e., macroscopic quantum 
tunneling. Some experimental studies of electric trans- 
port in superconducting wires at low temperatures have, 
indeed, demonstrated finite residual resistivity in suffi- 
ciently thin wires [2|, yj, consistent with the quantum 
hase slip (QPS) picture [J]. Other experimental studies 
however, did not confirm these results, suggesting, 
in particular, that the experiments in Ref. [2| could be 
explained in terms of thermally activated phase slips in 
the presence of Coulomb blockade. 

It has been argued in the literature [g, [7| that the 
controversy (of whether the observed phase sfips are of 
quantum origin) can be clarified by observing the coher- 
ent effects in fluctuations of electric currents, - the ther- 
mally activated phase slips, obviously, are not supposed 
to exhibit any phase coherence. In particular, based on 
phenomenological arguments, it was proposed that short 
and thin superconducting filaments can possibly operate 
as sources of coherent QPSs p-lsj, just like Josephson 
junctions in flux qubits. The experimental evidence for 
coherent QPS in Josephson junction arrays has recently 
been reported in Ref. [9]. Yet, no microscopic theory 
of coherent QPS-driven current fluctuations has been de- 
veloped so far. 

In this paper we provide a microscopic theory of co- 
herent QPS in a small superconducting (nano) ring. We 
argue that if a half-quantum of external magnetic flux 
is applied to the ring, it will exhibit coherent switching 
between states with different persistent currents, i.e., co- 
herent QPS transitions. The rate (frequency) of such 
QPS is controlled by the strength of the impurity poten- 
tial: For a translationally invariant system, i.e., in the 
"clean" limit, the QPS are suppressed due to the Landau 



criterion lfllll|. Indeed, a QPS event is accompanied by 
the production of sufficiently large momentum, which (at 
zero temperature) cannot be absorbed by the quasiparti- 
cles unless the current (or superfluid velocity) exceeds its 
critical value. In the presence of impurities the momen- 
tum conservation is "violated" and so the QPS can occur. 
However, as we show in this paper, the rate of the QPS 
is proportional to the Born's scattering cross-section of 
the condensate at the impurities, which is small for large 
QPS momenta, i.e., for wires whose diameter significantly 
exceeds the interparticle distance. Moreover, the value of 
the QPS frequency turns out to be strongly sample de- 
pendent: the ratio of its root-mean-square deviations to 
its average value is a constant which does not decrease 
with system's size, as in the case of incoherent QPS in 
sufficiently long wires 10|, [l^ . 



Before proceeding to the evaluation of the QPS fre- 
quency for a BCS superconductor, let us evaluate this 
quantity in the framework of the Bose-Luttinger liquid. 
The calculation for this model is rather straightforward 
and serves as a reference point for the more complex su- 
perconducting case. The Hamiltonian for such quantum 
fiquid confined to a ID ring with finite circumference L 
can be written as iJ = Hq + Hi UM , where 



Hf) — Hr 



Hj = 



]ujkalak + iEj/2){J-Jof (1) 



kjiO 



is the Hamiltonian of the liquid in the absence of impu- 
rity potential, i.e., for a translationally invariant system, 
and Hi is due to the impurity potential. In Eq. ([T]) cok 
is energy of a quasiparticle with momentum 2'Kk/L and 
al (ofc) are bosonic operators that create (annihilate) the 
quasiparticles. The Hj term in Eq. ([T]) is the energy as- 
sociated with the superfluid motion of the liquid and is 
characterized in terms of winding numbers J = 0, ±1, ..., 
i.e., eigenvalues of operator J. Constant Ej is related to 
the effective inductance of the liquid (typically of kinetic 
origin for both neutral systems and ultrathin supercon- 
ducting wires) and Jo is the "external" winding number 
supplied by the magnetic field (for charged particles) or 



by the rotation (for neutral particles). Note that in Eq. 
([1]) the states with different J's are decoupled from each 
other owing to the Landau criterion. 

The presence of Hi = J^ dxV{x)n{x), where V{x) is 
impurity potential and n is ID particle density, intro- 
duces a matrix element between states with different J's. 
This can be seen by writing n in terms of the displace- 
ment phase Q [13[ according to the relation n = no + 5n, 



Sn ~ Y.'yJ'Y exp{iie) + {J+y exp(-i/e), (2) 
;>o 

e = 2nnox + ^ (^) '^' (4 -|- a^,)e'"'^^/\ 

where operators J raise/lower the winding numbers J 
by ±1 , and a is the Tonks parameter, a = imo/{mc), m 
is particle's mass and c is sound velocity. In this paper we 
set both the Plank constant and the speed of light equal 
to 1. The value of (0|ifi|l), where |0) and |1) are the 
states of the liquid with zero and one winding numbers 
and zero quasiparticles, can now be easily evaluated with 
the use of Eq. ([2]). Retaining I — 1 terms only in Eq. 
([2]), we obtain 



\Hi\l) 




dxV{x)e 



27TinQX 



■ Z^fc 



(3) 



If the value of this matrix element is small compared 
to the gap between different quasiparticle states Aw = 
2ttc/L, the system's dynamics at low enough temper- 
atures (i.e., for T <^ Aw) reduces to the transitions 
between J = and J = 1 "ground" states (here 
and in the following we are interested in the degener- 
ate Jo = 1/2 case, which can always be achieved by 
the appropriate choice of magnetic field or rotation fre- 
quency). That is, if the system is initially set in a state 
with, say, J = 0, its probability to occupy that state, 
|(0(i)|0(0))p, will exhibit coherent oscillations with fre- 
quency Wo — |(0|iJi|l)|, as in the case of a two-state sys- 
tem. Assuming the Bogolubov's dispersion relation Wfc = 
c[(27rfc/L)2 + (27rfc/L)4CV4]i/^ where ^ = {mc)-'^ is co- 
herence length, we obtain that J2n ^n^ — ^'^ {'^^^^'^L/ttS,) 
in Eq. (O and 



^o^|Ko|x(^e/2ei/2i)^ 



(4) 



where Vn^ — J dxV{x)e'^'^^""^. It is interesting 
that Wo does not self-average with the growth of 
L. Indeed, if we assume that the disorder poten- 
tial is short range and Gaussian with {V{x)V{x')) = 



Vo'S{x 



we see that the average frequency (wo) = 



(7ri/V2)noVb2.^/2(7r^/2ei/2L)", while its root-mean- 
square deviations ((Awq))^/^ ~ (wq). Therefore we pre- 
dict that the frequency of coherent oscillations is strongly 
sample-dependent. While this fact may seem surprising, 
it is a direct consequence of systems 's coherence: Indeed, 



for a sufficiently large system, e.g., Refs. lOJjbJl, phase 
slips occur at different, disconnected from each other re- 
gions, which effectively leads to self-averaging. In our 
case the phase slips are strongly overlapping (note that 
Wo ^ c/L) and thus averaging over the disorder does not 
occur. 

Let us now consider the fermionic model. The Hamil- 
tonian density for electrons within the standard BCS 
model can be written as 



+e^ 



H = :^^i(*V - eAo)2V', 4 
Zm 

^3^ V4(r)V>^(r)^t^(rOV>,,(rO 



y(r)^t^, (5) 

f (AV'^V^i + h.c), 



where repeated index a means summation over f and J, 
spin components of fermionic field ipa- (r) . In Eq. ([5]) Ao 
is vector potential due to an external magnetic field and 
V{r) is impurity potential, to be specified below. The 
electrons are restricted to move in a wire of cross section 
Sq = Trrg (ro is assumed to be much smaller than the su- 
perconducting coherence length ^s), which has the shape 
of a ring (or torus) with circumference length L ^ ^g. 
We assume that the magnitude of the superconducting 
order parameter A(r) = |Ao|e'^^'''' is fixed and allow for 
variations of its "soft" phase 9 only. Such assumption 
is justified by the existence of the gapless (i.e., sound- 
like) mode arising in a quasi-lD situation (see below). 
In Eq. ([5]) we neglect the geometric inductance of the 
loop, which is small compared to the kinetic inductance 
for ultrathin wires. 

Unfortunately it does not seem possible to directly in- 
troduce the displacement phase O for the Cooper pairs 
as it was done in the previous model. For the conven- 
tional phase 9, however, one can derive an effective low 
energy action within the perturbative approach, i.e., as- 
suming that gradients of 9 as well as its time derivatives 
are small. In order to do so we apply a gauge transfor- 
mation -00- — >■ -00-6*^/^ to the Hamiltonian in Eq. ([5]) and 
expand the corresponding action in terms of 5t-0 and V9. 

The first order contribution to the effective action for 
the phase 9 is the Berry phase term [1J| 



Si^i dTd^r{i^l{r)i;^{r))drt 



(6) 



In Eq. ^ the averaging is taken with respect to the 
ground state of the electrons, but not with respect to the 
disorder potential V. The importance of par ticle density 
variations has been discussed in Refs. [id . Ill| . where 
it was shown that the QPS processes are suppressed by 
the Berry phase, e.g., Eq. ([S]), in translationally invari- 
ant systems. Such suppression is, essentially, a conse- 
quence of Landau criterion, discussed above. Note that 
the quantity J dTdr9 is a topological invariant: It is de- 
termined by a discrete set of points corresponding to the 
positions of the QPS rather than by the particular de- 
pendence of on T and r; see below. 



The second order terms in drO and V9 have been de- 
rived in Refs. |J|. For consistency we briefly outhne the 
calculation. First one introduces the scalar potential (j) by 
decoupling the Coulomb interaction term in Eq. ([S]) via 
Hubbard-Stratonovich transformation. The second order 
terms obtained after the averaging over the tp fields are 

82= I drSr [iy{idrO + 2e(j)f (7) 

+ {ps/2m){ye + 2eAo)' + {V(j)f /%^] , 

where z/ is free electron density of states and ps is the 
density of superconducting electrons. Note that Eq. ([7]) 
is also averaged with respect to the impurity potential. 
Such averaging is valid since is a slowly varying quan- 
tity. Moreover, since the thickness of the wire is small 
compared to the superconducting coherence length, is 
effectively a ID field. The field cj) can now be excluded 
from Eq. ([7]). In doing so one should keep in mind that 
(j) is also nonzero in the region outside the wire. The con- 
tribution of the outside region can be shown to dominate 
over the inside region and the last term in Eq. ([7]) can be 
written as J dxC(j)^ /2, where C ~ e[21n(L/ro)]^"'^ is the 
capacitance (per unit length) of the ring, e is dielectric 
constant of the surrounding medium and the scalar po- 
tential (j) inside the wire is assumed to be constant along 
the wire's cross section. Then ip can be integrated out 
and from Eqs. ([HI [7]) we obtain an effective Lagrangian 
density for 9{t,x): 



C 
£ = in(x)drO -f -—^{dr 
2e^ 



PsSg 

2m 



{dj + 2eAo) 



(8) 



In Eq. ^ we have defined the ID density n{x) — 
J d^r^p{r), where p{r) — {^pl{r)^p„{r)), and assumed 
that c^vSq ;» C. Aq in the last term in Eq. ^ is 
the component of the external vector potential Ao along 
the direction of the wire. The last two terms in Eq. ^ 
describe the Mooij-Schon gapless plasmon mode propa- 
gating with phase velocity Cs — (ps-Soe^/mC)^/^ |^^4|. 

We are interested in the probability amplitude for the 
system to remain in state |0), i.e., with zero current. If a 
half of magnetic flux, $0 = / dSBo — § d:x.AQ, is applied, 
i.e., $0 -^ 7r/2e, the zero current state is degenerate with 
the state with one winding number, |1), i.e., with electric 
current equal to 2TTeps/mL. The amplitude (0(t)|0(0)) 
can be evaluated using instanton method. The leading 
contribution is 



(0(t)|0(0)) r^l + K dn 

Jo Jo 



dro 



(9) 



X 



Jo Jo 



where action Sd is evaluated along the classical trajec- 
tory, which starts in state |0), then, at time ti passes 
on to the state [1), and at time T2 returns to |0). The 



coefficient K is related to the determinant for the fluctu- 
ations around the classical trajectory. It is independent 
of the first, top olog ical term in Eq. ([5]), and was esti- 
mated in Refs. [ll[: K ^ a^/(^^TQ), where tq = ^s/cs 
and as = 7r(ps5'oC/e^m)^/^. 

The classical trajectory action in Eq. ^ corresponds 
to a superposition of two phase slips (i.e., vortex at ti, xi 
and anti- vortex a.t T2,X2) [J]. For a finite ring we should 
account for the boundary conditions: The system is peri- 
odic in X direction (with period L) and infinite along the 
T direction (again we consider a zero temperature case). 
The (anti)vortex solutions for the phase 9 satisfying such 
boundary conditions has the form 



X tan 



X,t) = ±Sgn(T-Ti(2)) 

-1 



(10) 



g-^2,rc,|r-ri(2)l/i sin[27r(a; - a;i(2))/L] 



1 



,-27rCs|r-ri(2)|/-L 



<2)l/-^cos[27r(a;-Xi(2))/i^] 



Note that in the limit L -^ 00 we recover the well known 
expression 0{x,t) -^ ±tan~"'"[(a; — xi)/cs{t — ri)], e. g., 
Refs. Q. Using this solution and Eq. ^ we evaluate Sd 
in Eq. ^. After some calculation we find 

/•xi 

Sci = 2TTi dxn{x)+T]A<i>AT + 2as\n{L/S_s) (H) 

J Xo 



-2a, In 1 



2e 1- cos 



where At = T2 — ti, Aa; = X2 — xi, rj ^ AiTepsSo/mL 
and A* = $0 - 7r/2e. The first term in Eq. ([TT|) is due 
to the Berry phase in Eqs. ([H H]); it is controlled by the 
relative position of the vortices, but not by their internal 
structure. The second term in Eq. ()lip is proportional to 
the bias energy ?7A<I>, i.e., the energy difference between 
states with and without current. The last two terms are 
the the vortex self-energy and interaction energy respec- 
tively. 

At the point of degeneracy (for A<I> = 0), the inter- 
action energy between vortices vanishes when At = 00. 
Then the last term in the rhs of Eq. © is equal to 
a;gT^/2, where 



^l/2g-a,ln(L/C,)| / ^^f ^2^^ ff dx n{x) ^ ^ (^3) 



CJQ 



Moreover, it is straightforward to verify that the next, 
four-vortex contribution in Eq. ([5]) is ^ (wot)'*/4!, etc., 
and therefore (0(t)|0(0)) — cos{u!Qt), i.e., the system ex- 
hibits coherent oscillations with frequency ujq. 

For a translationally invariant system (with V — 0), 
n(x) — const and therefore the integral over x' in Eq. 
P^ vanishes and ujq = 0, as expected. For not too strong 
disorder and not too thin wires we may write n{x) = 
uq + Sn{x), with uq ^ 6n. Then 



dx 



I 27ri noa:;'+rf dx6n{x)\ 



1 



27rno 



d^ve 



2'KinQX 



Spir\ 



where p(r) — pQ + Sp{r). Since Sp{r) ~ / d^r'x(r — 
r')y(r'), where x is the static electronic susceptibihty, 
Eq. (fT^ is basically equivalent to that for the bosonic 
model, e.g., Eqs. ([31 0]). Note that the vortex self- 
interaction term — asln(L/^s) in Eqs. pTlfT^ coincides 
with the last exponent {^ J2n^n^) ii^ ^1- ®j arising 
from the evaluation of (exp (iQ)) in the bosonic problem. 
Again, as in the bosonic problem, we expect that 
the QPS frequency is strongly sample-dependent, i.e., 
it is a random number whose average is comparable to 
its root-mean-square deviations. In order to estimate 
the typical (average) frequency of QPS oscillations, we 
must evaluate {\ J dxe^""^"^ Sn{x)\'^). We assume that 
the disorder is due to random impurities carrying charge 
e, i.e., V{r) = J2i=i^^/\^ ~ ''^il- The response func- 
tion X for the Hamiltonian in Eq. ([5]) is generally not 
well known at finite wavevectors. However, the momen- 
tum 2TmQ = 2npoSo transferred from the condensate 
to the impurities as a result of a QPS event is rather 
high compared to both the Fermi wavevector kp and 
the inverse screening radius. At such high momenta 
the electronic response function is that of free electrons, 
x{q) — Xo~*^°°('?) — — 4mfc|,/(37r^(7^), and therefore 



^3j.g2^inox 



Spir)f) = [x^i27rno)eyTrnlfm, (13) 



where Ni is the total number of impurities in the sample. 
The right hand side of Eq. (fT3|) is ^ Sq and therefore 
the QPS rates are strongly suppressed in wires, whose 
diameter significantly exceeds the Fermi wavelength. 
Then we obtain an estimate for the average frequency: 



(^o) 



2 Arl/2 



(14) 



where as is Bohr radius. Let us estimate the value of 
(wo) for typical experimental parameters. We consider a 
MoGe wire, e.g. Ref. [2|, of diameter 2ro = 5nm and 
length L = 1 fj,m. The density of conduction electrons 
po can be estimated from the data on conductivity in 
the normal state, e.g. Ref. [15[, a = e'^po/kpl- For an 



amorphous material, such as MoGe, the mean free path 
is of the order of interatomic distance / ~ 4A, e.g. Ref. 
[15| , and thus we obtain po — 10^* m~'^. In order to es- 
timate as and Cg we need values of ps and C. Assuming 
that the wire is covered by a dielectric material, e.g. Ref. 
0, S Qj we estimate that C ^ I. The density of su- 
perconducting electrons is related to the normal electron 
density asps — Poi^UsY \^- Taking ^s — Snm, e.g. 
Refs. [3, ll5|, we obtain ps ~ lO^'^m"'^ and a^ — 0.5 and 
Cs ~ 7 X 10^ m/s. Finally, the number of impurities Ni is 
of the order of the total number of atoms in the wire and 
so we take iV^ ~ 3 x 10^. With such parameters we ob- 
tain (wq) ^ 10^ s~^. As we pointed out above, the QPS 
rate is significantly suppressed by the small "scattering 
amplitude", e.g. Eq. p^ . due to the high value of the 



QPS momentum. The value of (wg) can greatly increase 
in thinner wires, e.g. Eq. (J14p : For wires with diame- 
ter 2nm we obtain (uq) ^ 10^ s^^. The coherent QPS 
in such wires can presumably be detected by inductively 
coupling the ring to a SQUID. 

In summary we have studied a possibility for gener- 
ation of coherent QPS in small superconducting rings. 
We have evaluated the frequency of such transitions and 
found that it is appreciable only in ultrathin wires, i.e., 
with diameters not too strongly exceeding the electron's 
Fermi wavelength. We have also shown that such fre- 
quency is a non-self-averaging, sample-dependent quan- 
tity. 
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